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ABSTRACT 

We measure the location and energetics of a S iv BALQSO outflow. This ouflow has 
a velocity of 10,800 km s“^ and a kinetic luminosity of 10^® ’’erg s“^, which is 5.2% 
of the Eddington luminosity of the quasar. From collisional excitation models of the 
observed Siv/Siv* absorption troughs, we measure a hydrogen number density of 
nH=10^'^ cm“^, which allows us to determine that the outflow is located 110 pc from 
the quasar. Since S iv is formed in the same ionization phase as C iv, our results can 
be generalized to the ubiquitous C iv BALs. Our accumulated distance measurements 
suggest that observed BAL outflows are located much farther away from the central 
source than is generally assumed (0.01-0.1 pc). 

Key words: galaxies: quasars - galaxies: individual (SDSS J083126.15-1-035408.0) - 
line: formation - quasars: absorption lines 


1 INTRODUCTION 


Broad absorption line (BAL) outflows are detected as 
absorption troughs that are blueshifte d in the rest- 
fram e spectrum of 20-40% of qua s ars llHewett fc Folt j 


20031: Ganguly fc BrothertonI l2008l : iKniAe et al.l ^0^ 


Dai. Shankar fc Sivakofll l2008h . From their detection rate, 
we deduce that these outflows cover on average ~20-40% 
of the solid angle around the quasar. Such large opening 
angles allow for efficient interaction with the surrounding 
medium. As shown by simulations, the mass, momen¬ 
tum and especially the energy carried by these outflows 
can play an important role in the evolution of galaxies 
and their environments fe.g. |ScantiaBigCQ_^Llh| [loo^; 
Levine&_Gnedi^200^ Ho 2 ldnset_al 2006: Cattaneo_^£_al 


20091 : Ciotti. Ostriker fc Progal l2009l. l2010l: 

20101 : 


Gilkis fc Sokeil 20121 : Choi et ^ 2014 1. Theoretical 


Ostriker et al 


studies show that such interactions can provide an ex¬ 
planation for a variety of observations: the self-regulation 
of the growth of the supermassive black hole and of the 
galactic bulge, curtailing the size of massive galaxies, 
and the chemical enrichment of the intergalactic medium. 
Thes e processes are pa rt of the so-called AGN feedback 
(e.g. ISilk fc Reed 199^: Di Matteo. S oringel fc Hernauis^ 


20051: 


jGnrmaimBarai fc hlart^ 


Hopkins. Murray fc ThompsonI 20091 : lElvij 

Zubovas fc Navakshin 20141 . and references therein). 


20091: 


20061 : 


The importance of quasar outflows to AGN feedback 
depends on the mass-flow rate (M) and kinetic luminosity 
(E]z) of the outflowing material. An value of at least 0.5% 


Table 1. Physical properties of energetic quasar outflows. 


Object 

V 

(km/s) 

R 

(pc) 

M 

(Mo/yr) 

Log Ek 
(erg/s) 

-^fc/Tgdd 

(%) 

SDSS 40831+0354^ 

-10800 

110“ 

135 

45.7 

5.2 

SDSS JIIO6+I9392 

-8250 

320“ 

390 

46.0 

12 

HE 0238-19043 

-5000 

34005 

140 

45.0 

0.7 

SDSS 40838+2955^ 

-5000 

3300“ 

300 

45.4 

2.3 

SDSS 40318-06005 

-4200 

6000“ 

120 

44.8 

0.13 

R from: “high-ionization Siv*/Siv: 

; “high-ionization Olv*/OlV: 


^low-ionization Siii*/Siii 

^^See Section l5.1l for I/Edd determination. _ _ 

Referenc es. ( 1) this work; (2) iBorguet et al.l |201^: (3) 

][ 2 ^; 


lArav et al. I I2013I: (4) iMoe et 151 iDunn et'al]|20ld . 

dHopkins fc Elvisll2010l l or 5% dScannapieco fc (5hll2004h of 
the Eddington luminosity (I/Edd) is deemed sufficient to pro¬ 
duce the aforementioned feedback effects. A crucial param¬ 
eter needed to determine Ek is the distance R to the out¬ 
flow from the central source. Lacking spatial image infor¬ 
mation, we deduce R from the value of the ionization pa¬ 
rameter {Un, see Section[4|) once the hydrogen number den¬ 
sity (wh) of the gas is known. Our group has determined 


(e.ff.lMoe et al.ll2009l: Dunn et al.ll201C 

:lBautista et al.ll2010l: 

lAoki et al.l 201ll: Boreuet et al.l 2012b 

1 by utilizing absorp- 


tion lines from excited states of singly-ionized species (e.g. 
Fell* and Sill*). 

Absorption troughs from singly-ionized species classify 
an outflow as a LoBAL. The lower detection rate of Lo- 
BALQSO in spectroscopic surveys (3-7%) compared with 
20-40% for Giv BALQSO dDai. Shankar fc SivAofj|2012f l 
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raises the question of whether the determinations obtained 
for these objects are represen tative of the ubiqu itous high- 
ionization Civ BALQSO fsee lPunn et al.ll201(ll l. The most 
straightforward way to avoid such uncertainty is to observe 
outflows that show absorption lines from excited states of 
ions with a similar ionization potential to that of C iv. The 
optimal io n for ground-based observations is S iv (see dis¬ 
cussion in iDunn et al.ll20l3 : lArav et"^ ][2^. 

To realize this, we conducted a survey using VLT/X- 
shooter between 2012 and 2014 aimed at finding quasar out¬ 
flows showing absorption from Siv A1062.66 and the excited 
state SIV* A1072.97. Of the 24 objects observed, two ob¬ 
jects (SPSS J1106-I-1939 and SPSS J1512-I-1119) have been 
published (iBorguet et al.l l2013l . hereafter Paper I). SPSS 
J1106-1-1939 yielded the most energetic SiV BAL outflow 
to date with erg s“^ (see Section 15.31 for fur¬ 

ther discussion), whereas the SPSS J1512+1119 outflow has 
10'^®'® erg s“^. In this paper, we present the analysis of 
an outflow from SPSS J0831-I-0354, which exhibits similar 
properties to that of SPSS J1106-I-1939. This includes the 
presence of P v and Siv/Siv* troughs, thereby allowing us 
to determine the distance and energetics (see Section [5. Ill of 
the SPSS J0831-I-0354 outflow. Table [1] summarizes the cur¬ 
rent state of the held by listing the energetics of prominent 
outflows analyzed by our group. 

The plan of this paper is as follows. In Section [2] 
we present the VLT/X-shooter observations of SPSS 
J0831-I-0354. In Sectlon[3]we identify the absorption troughs 
of the outflow and measure the ionic column densities which 
are used to determine the number density (Section 13.411 and 
the photoionization solution (Section|3)). In Section[S]we de¬ 
termine the outflow distance and energetics and compare 
them to outflows from other objects. We summarize our 
method and findings in Section |5] 


2 OBSERVATIONS AND DATA REDUCTION 

SPSS J0831-f0354 (J2000: RA=08 31 26.15, PEC=+03 54 
08.0, z=2.0761) was observed with VLT/X-shooter {R ~ 
6000-9000, see Paper I for instrument specifications) as part 
of our program 92.B-0267 (PI: Benn) in January 2014 with 
a total integration time of 10640s. 

We reduced the SPSS J0831-I-0354 spectra in a similar 
fashion to those of SPSS J1106-I-1939 (detailed in Paper I): 
we rectified and wavelength calibrated the t wo-dimensional 
spect ra using the ESO Reflex workflow dBallester et al.l 
[2012), then extracted one-dimensional spectra using an op¬ 
timal extraction algorithm and finally flux calibrated the 
resulting data with the spectroscopic observations of a stan¬ 
dard star observed the same day as the quasar. The one¬ 
dimensional spectra were then coadded after manually per¬ 
forming cosmic-ray rejection on each spectra. We present 
the reduced UVB-I-VIS spectrum of SPSS J0831-I-0354 in 
Figure [T] 


3 SPECTRAL FITTING 

Absorption troughs associated with Hi, Civ, Nv, Mgii, 
Aim, Sim, Siiv, Pv and Siv/Siv* ionic species are seen in 
the spectrum (see Figure [TJ. The high velocity (u=—10800 


km s“^) and width (3500 km s“^) of the C iv trough satisfies 
the definition of a BAL outfl ow (i.e. Umax > 5000 km s~^ 
and width > 2000 km s“^, see lWevmann. Carswell fc SmithI 

a . The balnicity index of the C iv trough is 1100 km s“^ 
N 


TWevmann et al.lll99lh . 


3.1 Unabsorbed Emission Model 


We co rrect the spectrum for galactic extin ction {E{B — V) — 
0.025: ISchlegeh Finkbeiner fc Pavi^ll998ll using the redden¬ 
ing curve of Cardelli. Clayton fc Mathij lll98^ . We then 
fit the continuum with a cubic spline resembling a power 
law of the form F{X) = F’iioo(A/1100)“, where Fiioo = 
1.58 X 10“^®erg s“^cm“^A“^ is the observed flux at llOOA 
(rest-frame) and a ~ —1.16. We model the broad emission 
lines (BEL) with a sum of one to three Gaussians, which 
does not significantly affect the column density extraction 
for most ions since the high velocity of the outflow (u=- 
10800 km s“^) shifts the BALs far from the wings of their 
corresponding BELs. However, the S iv and S iv* BALs lie 
within the Ovi BEL, which we model by scaling the Civ 
BEL template to match the peak emission of the Siv BAL 
region. 


3.2 The Blended Troughs 

We model the optical depth of the absorption troughs with 
three Gaussians 

3 

t{v) = ^ Ti exp 

with the same centroid and width in all the observed 
ionic troughs. Our first Gaussian (t)i=—12,100 km s“^, 
FWHM=1460 km s“^) is chosen specifically to match the 
blue wing of the Siiv A1393.75 line, as well as the blue 
wing of the Aim A1854.72 line. The second Gaussian 
(u 2 =—10,800 km s“^, FWHM=940 km s“^), which rep¬ 
resents the majority of the absorption, targets the Siiv 
and Pv doublets. A third Gaussian (t) 3 =—10,000 km s“^, 
FWHM=590 km s~^) is needed to match the red wing of 
the Aim A1862 line, as well as the red wing of the C iv and 
N V blends. These three Gaussians form the template used 
to extract the column density (following the same procedure 
as Paper I) from the observed troughs of all ions except for 
Pv and Sim, which will be discussed in the next section. 
For the three-Gaussian template, we use a covering factor 
(Cion) that is constant throughout the trough (i.e. for all 
three Gaussians), and varies only between ions . For a de¬ 
script ion of the partial covering (PC) model, see lArav et al.l 
ll2008ll . For ions detected by a singlet, we take the apparent 
optical depth (AOD) measurement as a lower limit. 


FWHM = 2V21n2a (1) 


3.3 Velocity-dependent covering: Pv and Sim 

Although the central Gaussian of our template (see Sec¬ 
tion 13.211 can be used to model the P v and Si iv doublet 
troughs, the ratio between the red and the blue compo¬ 
nents of each doublet approaches 1:1. This suggests that the 
troughs are strongly influenced by the velo city-dependent 
partial covering of the emission source (see iBorguet et al.l 
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Figure 1. VLT/X-shooter spectrum of the quasar SDSS J0831+0354 (z=2.0761). We label the ionic absorption troughs associated with 
the outflow, and represent the unabsorbed emission model with the red dashed line (see Section l3.1ll . Narrow absorption from intervening 
systems appear throughout the spectrum, and terrestrial absorption from molecular O 2 in our atmosphere is seen near 6850A and 7600A 
observed frame, but none of these features affect the analysis presented here. 


l2012al : lArav et al.l 11999^ 1. Due to the presence of Lya for- 
rest intervening troughs, we proceed by modeling the nor¬ 
malized flux of the doublets with a smooth function. For 
phenomenological flexibility, we built a function from the 
product of two logistic functions (hereafter logistic fit) 


1 ' (1 _|_ glUrCv-Vr)) ^ ■' 

where A is the maximum depth of the logistic fit, vt, and v-r 
are the velocities of the blue and red wing half-maximums 
respectively {vb < Ur), and Wb and Wr are the slope/variance 


of the blue and red wings respectively {wb < 0 < Wr)- This 
function is particularly suited to fitting asymmetric troughs. 

We model each P v doublet component (normalized and 
in velocity-space) using a logistic fit, and determine the 
velocity-dependent covering fraction and opti cal depth so¬ 
lution (using Equation (3) of iDunn et al.ll2010l b Integrating 
this optical depth over the trough yields the Pv column 
density for the partial covering model given in Table (2] 

Sim is a singlet, thus we cannot determine a covering 
model for this ion. Therefore, we use the covering factor of 
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Figure 2. Fits to the absorption troughs observed in the X-shooter spectrum of SDSS J0831-I-0354. Each trough is fit by scaling the 
optical depths of three template Gaussians whose centroid and width are fixed among the ions (see Section l3.2l l. The Gaussian templates 
for each transition are shown in the same color (red or blue), and the combined absorption from blended lines is represented as a black 
solid line. The first row of plots show the apparent optical depth (AOD) fits to Sim (and Hi), Mg II and Aim. Due to the influence of 
non-black saturation in the Si IV and Pv troughs (see text), we show the velocity-dependent partial covering (PG(v)) fit to those ions. 
The plot in the lower-right shows the PC(v) fit to the S IV and S IV* troughs in velocity space. The fits are spanned by shaded contours 
representing the upper and lower errors assigned to the fits. 


Si IV as a proxy to extract the Si ill PC column density given 
in TableO since Si iv is the next ionization stage of the same 
element. The lower error of 0.11 dex for the adopted Sim 
measurement includes systematic errors in continuum place¬ 
ment and blending of the intervening Lya forest. We note 
that, by coincidence, the difference between the PC measure¬ 
ment and the apparent optical depth (AOD) measurement 
is also 0.11 dex. 

3.4 The density-sensitive tronghs: S iv and S iv* 

The column densities from S iv and its excited state S iv* 
play an essential role in our analysis: the ratio of level pop¬ 
ulations between S iv* and S iv is dictated by collisional ex¬ 
citation, which depends on rie and is relatively insensitive to 
temperature (see Figure 8 of Paper I). 

As in the case of Sim, Siv A1062.66 and Siv* A1072.97 
are both singlets (the lines arise from different energy levels 


with a priori undetermined populations), thus a velocity de¬ 
pendent covering solution cannot be determined from these 
two lines alone. Therefore, we follow the same method as 
with Sim, using the velocity-dependent covering template of 
Si IV to model the S iv/S iv* troughs with the PC(v) method. 

We begin this method by modeling the Siv/Siv* 
troughs using a logistic curve that closely resembles the cen¬ 
tral Gaussian of the three Gaussian template. This resemb- 
lence ensures that the S iv absorption model is kinematically 
related to the other ionic troughs in the outflow. The logistic 
curve was fit to the S iv troughs by scaling the amplitude un¬ 
til it fit the broad absorption features in the spectrum while 
excluding the narrow Lya troughs. To assign the upper and 
lower errors to this fit we again scale the depth of the lo¬ 
gistic template until the resulting template clearly over- or 
under-predicts significant portions of the broad absorption 
trough. The quality of this fit and of the stated errors can 
be judged in the lower right panel of Figure [2] 
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Figure 3. Density diagnostic using the Siv*/Siv ratio for the 
SDSS J0831+0354 outflow, presented in the same manner as Pa¬ 
per I: we plot the S IV*/S IV theoretical ratio versus electron num¬ 
ber density and overlay our measured value (for two absorber 
models, AOD and velocity-dependent PC(v)) to determine the 
number density and its errors. 

The optical depth of the fit (and its errors) was trans¬ 
lated to Siv/Siv* column density measurements using the 
same procedure as in Sim (i.e. adopting the velocity- 
dependent covering solution for Siiv). These column densi¬ 
ties are reported in the fourth column of Table [2] A similar 
value (differing by only 0.05 dex) was obtained with the P v 
velocity-dependent covering template, demonstrating the ro¬ 
bustness of this method. 

To obtain the AOD column density, we follow the same 
procedure as the PC(v) modeling, instead using the central 
Gaussian of the three Gaussian template described in Sec¬ 
tion For completion, we also apply the constant partial 
covering template of Silv to the Siv/Siv* absorption (us¬ 
ing the central Gaussian as the optical depth template), and 
report the PC column density in the third column of Table[2] 

We deduce the rie of the outflow by comparing the 
Siv*/Siv column density ratio to predictions (s ee Figure [3l) 
made with the Chianti 7.1.3 atomic database (iLandi et al.l 
l2013h . Using the column densities for S iv and S iv* reported 
in Tableland the electron temperature found for our best- 
fit Cloudy models from Section 2] (9000 K, a weighted av¬ 
erage for Siv across the slab), we find log(ne) = 4.42(to'22 
cm“® by averaging the two absorption models (reported in 
Table 2]) and adding their errors in quadrature. 


4 PHOTOIONIZATION ANALYSIS 


We use photoionization models in order to determine the 
ionization equilibrium of the outflow, its total hydrogen col¬ 
umn density (Ah), and to constrain its metallicity. The ion¬ 
ization parameter 


_ Qh 
dvrA^criH ’ 


( 3 ) 


Table 2. SDSS J0831-I-0354 log column densities (cm“^) tor 
the three absorption models AOD, PC and PC(v) (velocity- 
dependent covering). Upper limits are shown in red and lower 
limits are in blue. 


Ion 

AOD 

PC 

PC(v) 

Adopted“ 

Hi 

>14.96 

>15.17 

— 

>15.17 

He I* 

<15.14 

-- 

— 

<15.14 

Cii 

<14.46 

— 

— 

<14.46 

Civ 

>15.50 

>15.65 

— 

>15.65 

Nv 

>15.60 

>15.79 

-„ 

>15.79 

Mg II 

>13.78 

13.87+°« 

— 

13.87+°^ 

Alii 

<13.07 

— 

— 

<13.07 

Aim 

>14.10 

14.27«13 

— 

14 

13 

Sim 

>14.27 

>14.35 

>14.38 

>14.38^-0.11 

Si IV 

>15.17 

>15.44 

>15.40 

>15.40 

Pv 

>15.18 

1 c oq+0.05 

iO.OO_Q Q5 

15 47+0-15 

10.^1 _Q J^Q 

15 47+0-15 

Siv 

>15.53 

15.48«:?^ 


15-60+°:r5 

Siv* 

>15.35 

15.40«:?^ 

15.39t°:“ 

15.39l°;l° 


“Adopted value for photoionization modelling (see text). 

*’We use the difference between the Si III AOD and PC measure¬ 
ments as a lower error on the Sim lower limit. 


Table 3. Number density measurements. The values correspond 
to the horizontal error-bars in Figure [3] 



AOD 

PC(v) 

Adopted 

log(ne) 

4 47+*^-^^ 

4 

^•^^-0.22 


(where Qn is the source emission rate of hydrogen ionizing 
photons, R is the distance to the absorber from the source, c 
is the speed of light, and nu is the hydrogen number density) 
and Ah of the outflow are determined by self-consistently 
solving the ionization and thermal balance equations with 
version cO S.Ol of the spe c tral s ynthesis code Cloudy, last de¬ 
scribed in iFerland et al.l (l2013l l. We assume a plane-parallel 
geometry for a gas of constant tih and initially choose solar 
abundances and the UV-Soft SED which is a good repre- 
sentation of radio -quiet quasars (described in Section 4.2 of 
iDunn et al.ll20l3 . other SEDs and metallicities will be ex¬ 
plored in Section 14.21 here!. Eor the chosen SED and metal¬ 
licity, we generate a grid of models by varying Ah and Uh- 
Ionic column densities (Aion) predicted by the models are 
tabulated and compared with the measured values in order 
to determine the models that best reproduce the measured 
Aion. For a more ela b orate descr i ption of this method, see 
iBorguet et al.l ll2013ll ; lArav et ^ (I2OI3I') . 


4.1 Photoionization Solution 

In Figure U] contours where model predictions match the 
measured Aion are plotted in the Ah — Uh plane. Note that 
Figure |4] presents models for two different metallicity val¬ 
ues (IZ 0 and 4Zq, see discussion in Section (4.211 . The ionic 
column densities used in the photoionization modelling are 
listed in the last column of Table consisting of the PC 
measurements where available and with PG(v) prioritized 
over PG. In FigureU) we show only the ions which dominate 
the solution; the lower limits (Hi, Giv, Nv and Silv) and 
upper limits (Hel*, Cll and Alll) are trivially satisfied by 
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Log(UH) ° 


Figure 4. Phase plot showing the photoionization solution using 
the UV-Soft SED for gas with one and four times solar metallicity. 
Each colored contour represents (for one of the two metallicities) 
the locus of models (Uht^h) which predict a column density 
consistent with the observed column density for that ion. The 
bands which span the contours are the l-cr uncertainties in the 
measured observations. The dashed line indicates the Si III lower 
limit. For each metallicity, the black “X” is the best ionization 
solution and is surrounded by the contour (see text). 


the contraints of Al ill and S iv. The Mg ll contour lies be¬ 
tween the Aim and Siv contours, and is thus also satisfied. 

For Z = IZq, the Pv contour requires a higher Ah 
than the S iv contour in regions also satisfied by the Si ill 
lower limit (above the dashed line in Figure |4)). This over¬ 
prediction of the S IV column density will be discussed in 
Section [5. 2 1 The solution is also influenced by the Aim con¬ 
tour which lies parallel to S iv. This results in a poor fit 
(Xred=10.6), since no models can simultaneously predict the 
observed Siv and Aim column denisties. 


4.2 Dependence on SED and metallicity 

We consider the sensitivity of the photoionization solution to 
o ur choice of SED an d metallicity by following the approach 
of lArav et al.l ll2013t) . thereby allowing for comparison with 
previous outflows under the same set of assumptions. We 
find the ph otoionization solu tion in six different cases (three 
SEDs from lArav et al.ll2013l each with IZq and dZ©), plot 
the contours in the top panel of Figure [5] and report the 
solutions in Table O The change in metallicity from IZq to 
4Zq, scaled accord ing to the Cloudy starbur st schema (fol¬ 
lowing grid M5a of lHamann fc Ferlandlfl993l ) . has the same 
effect with each SED, decreasing the log(AH) of the solution 
by 0.4 dex. This is illustrated in Figure [T] which shows the 
photoionization solution for metallicities of both Z = IZq 
and Z = 4Zq. We will discuss alternative elemental abun¬ 
dances in Section 15721 here we consider the same metallicities 
used for other outflows studied by our group. 

The different SEDs we use spread the solution over 0.4 
dex in log((7H), which is comparable to the errors on each 
solution. In Table 2] we report Qh as well as the bolometric 
luminosity TboI for the three SEDs by fitting them to the 
measured flux (corrected for Galactic reddening) at 1100 A 
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Figure 5. (top) Phase plot showing the photoionization solution 
for three SEDs (HE0238, MF87 and UVsoft, see descriptions in 
lArav et al.ll20l3‘) and two metallicities (one solar (Zl) and four 
solar (Z4)) for a total of six models. Contours of equal -Ek/UEdd 
(assuming fixed tih) are shown as parallel, thin solid lines assum¬ 
ing the UV-Soft SED (see text), (bottom) Same as above, but 
with the horizontal axis converted to distance, which incorporates 
the errors on riu and Uu into the horizontal errors on R. 


Table 4. Properties of the chosen SEDs. 


SED 

HE0238 

MF87 

UV-Soft 

log(LBol) (erg s-l) 

46.79 

47.03 

46.87 

log(QH) (s“l) 

56.53 

56.95 

56.72 


(in the rest-frame) and integrating over the whole energy 
range. 

One important motivation to find both the ionization 
parameter C7 h and the number density tih is to determine 
the distance of the outflow R from Equation ®. We there¬ 
fore convert the ionization parameter log(t/H) to distance R 
using Equation and present the same six photoionization 
solutions in the bottom panel of Figure [5] The uncertainties 
on Uu and nu are incorporated into our value for R, which 
forms the horizontal spread of the contours in Figure[5] Since 
Ah is not involved in the distance determination, the verti¬ 
cal spread of the contours in Figure [5] remains unchanged. In 
both panels of Figure [51 we also show the contours of equal 
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Eii/L-Edd (see Section ISTTI) assuming the UV-Soft SED, not¬ 
ing that the same contours for the different SEDs are shifted 
in log(A^H) by no more than 0.1 dex. 


5 RESULTS AND DISCUSSION 
5.1 Energetics 

Assuming the outflow is in the form of a thin partial shell, 
its mas s flow rate (M) and ki netic luminosity (Ek) are given 
by fsee lBorguet et ai1l2012bl . for discussion) 

M = 4'KRQjj,mpNHV (4) 

Ek = 2'KRQ.^j.mpNHV^ (5) 

where R is the distance from the outflow to the central 
source, is the global covering fraction of the outflow, /r 
= 1.4 is the mean atomic mass per proton, rup is the mass 
of the proton, Nh is the total hydrogen column density of 
the absorber, and v is the radial velocity of the outflow. Us¬ 
ing the parameters reported in the preceeding sections, we 
calculate the energetics for the six models we have consid¬ 
ered for the S iv outflow of SDSS J0831+0354 and report 
the relevant values in the first six rows of Table [S] (the other 
four rows are comparison outflows that will be discussed in 
Section [5.31) . We adopt the UV-Soft SED with Z = 4Zq as 
our representitive model for this outflow, which is the same 
model chosen in Paper I, and yields a conservative mea¬ 
surement of Ek. As in Paper I, we use D = 0.08, which is 
appropriate for Siv BAL outflows. 

A s noted in the In t roduc tion, an Ek value of at least 
0.5% (iHopkins fc Elvid l2010l l or 5% dScannapieco fc OhI 
I 2 OO 4 I I of the Eddington luminosity (LEdd) is deemed suf¬ 
ficient to produce significant AGN feedback effects. We de¬ 
termine the Eddington luminosity Lsdd = 10^^'° erg s“^ 
from the ma ss of the SMBH (e.g. using Equation (6.21) of 
iKroIikl dlOOgl U. The mass of the SMBH is deter mined using 
the v irial mass estimator from Equation (3) of dPark et al.l 
I 2 OII given the caveat that the scaling relationship extrap¬ 
olated to high-z/luminosity quasars is not yet firmly es¬ 
tablished). This estimate requires the continuum luminos¬ 
ity at rest-frame 1350A and the FWHM of the C iv BEL, 
both of which are directly measured from our spectrum. 
We use this method to determine EEdd for the quasars 
SDSS J1106+1939 and SDSS J0838+2955. The Civ BEL 
in SDSS J0318-0600 is not well-defined, so we perform a 
similar analysis using the M gll BEL and Equation (1) of 
IVestergaard fc Osmei d2009l i. The only spectral coverage of 
HE 0238-1904 is in the extreme-UV, which does not cover 
any diagnostic lines sufficient for thi s method. We the refore 
use the assertion (from Section 9 of lArav et al.l[20l3 ') that 
the Eddington luminosity is approximately the bolometric 
luminosity for our comparison in Ta ble [B _ 

We note that recent studies dLuo et al.l 120141 ') have 
found X-ray weak SEDs for BALQSO that are much softer 
than the UVsoft SED we use here. To explore the effects 
of such an SED o n our results, w e constructed an SED 
from the object in ILuo et ^ d2014l ') showing the most ex¬ 
treme X-ray softness (PG 1254-1-047) by interpolating be- 
twee n the given flux points (see right-centre panel of Figure 
3 of iLuo et al.ll201-i L At 2 keV the flux is three orders of 
magnitude lower than in our UVsoft SED. The photoioniza¬ 


tion solution for this very soft SED leads to a distance of 
51 pc and Ek=10'*^'®erg s“^ for the outflow using Z—4 Zq. 
This distance and energy is lower by a factor of two and 
six respectively, compared to the values obtained with the 
UVsoft SED (110 pc, Ek=10^® ’^erg s'^) for Z^4 Zq. 

We caution that the cause of the observed X-ray weak¬ 
ness can be attributed to optically thick absorption from the 
outflow itself, rather than intrinsically weak X-ray emission 
or obscuration between the central source and the outflow. If 
this is the case, then the SED incident on the outflow would 
likely resemble one of the other SEDs presented here. It is 
the incident radiation that will determine the photoioniza¬ 
tion solution. 

5.2 The S IV discrepancy 

Our photoionization solution (for the UV-soft SED) over¬ 
predicts the S IV total column density by a factor of eight, 
which is not alleviated appreciably by changing the SED or 
metallicity. We offer two resolutions to this discrepancy: a 
two-phase solution, or individual abundance scaling of sul¬ 
phur. _ 

ll lArav et al.l ll2013l ~l showed that two ionization phases 
exist in some outflows (in that case smaller, denser cloudlets 
embedded in high-ionization gas with much larger Ah), and 
can account for the appearence of low-ionization species 
that are inconsistent with the ionization solution obtained 
from higher-ionization species. Siv and Pv are both high- 
ionization species, thus their column densities can be pro¬ 
duced in a separate phase than the low-ionization species 
(Mg II and Aim). The Siv and Pv contours in FigureUcon- 
verge at high ionization parameters (log [/h> 0.3). A high- 
ionization phase at this ionization parameter would satisfy 
both the Siv and Pv constraints, but under-predict the 
Mg II and Aim column densities, which would need to be 
satisfied by a low-ionization solution. A two-phase solution 
could resolve the S iv discrepancy, but such a solution is not 
well-constrained with the measurements obtained from our 
spectra. However, any high-ionization solutions that adhere 
to the P V constraint would move towards the upper-right 
corner of Figure (4] Since Ek oc Ah/VUh, this increase in 
both Uh and Ah by 0.6 dex would result in an increase in 
Ek of 0.3 dex. 

2) The discrepancy between S iv and the other ions 
may be alleviated by considering a higher metallicity, as 
outflows are known to have supersolar m etallicities (e.g. 
iGabel. Arav fc Kimll2006l : lArav et al.l[200^ L In our Gloudy 
models the discrepancy persists with different metallici- 
tie s since the abunda n ce sc aling schema from grid M5a 
of iHamann fc Ferlandl lll993ll does not increase the rela¬ 
tive [S/P] or [S/Al] abundances. An alternative enrichment 
model that results in [S/P]=-0.9 and [S/Al]=-0.9 for a cer¬ 
tain metallicity would decrease the S iv column density pre¬ 
dicted by the solution from the P v and A1 m constraints 
by a factor of eight, thus eliminating the discrepancy that 
appeared in the solar abundances case. Considering the ab¬ 
sence of a complete AGN abundance scaling model, this so¬ 
lution is not as implausible as it seems; e.g. iBallero et al.l 
ll2008ll derives relative abundances such as [Si/C]=0.83 for 
E = 7.22Z0. 

If neither of these scenarios apply, then we conclude that 
the S IV column density is indeed a factor of eight higher than 
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Table 5. Physical properties of energetic quasar outflows (object names abbreviated from Table[T}. 


Object, SED 

log(LBol) 

V 

log(Uif) 

log (Nh) 

log(ne) 

R 

M 

iog(Ek) 

Ek/L^dd 


(erg s-l) 

(km s“l) 


(cm“2) 

(cm“®) 

(pc) 

(Mgyr-i) 

(erg s-l) 

(%) 

J0831 HE0238 Z1 

46.8 

-10800 



^■^^-0.22 

78lt^ 

9qn+'^’^^ 

Z9U_i3o 

45.91“;^ 

o + ii 

^-4.4 

J0831 MF87 Z1 

47.0 

-10800 


99 9+ U-^ 

^•^^-0.22 

180l®® 

3901®™ 

46-11^ 

10 + 19 

J0831 UVsoft Z1 

46.9 

-10800 

-0.26t°;t 

22.51“® 

^•^^-0.22 

110l®“ 

4ioig[; 

46.2l“;4 

14+18 
1^—7 7 

J0831 HE0238 Z4 

46.8 

-10800 

-0.16+0® 

21 9+^® 

4 42'*"®-^® 
^•^^-0.22 

70 + 33 
'^-24 

76+120 

'^-46 

45 4+^’^ 

9 

^•’^-1.6 

J0831 ME87 Z4 

47.0 

-10800 


21.61“;® 

^•^^-0.22 

200l^^ 

note™ 

45 6+^-“^ 

^o.u_o 4 

q q+5.8 
3 -o_2 3 

J0831 UVsoft Z4 

46.9 

-10800 

-0.25l°j 

22.o1“;4 


notll 

isol^f 

45.71“-® 

c 9+4.8 
^•^-2.9 

J1106 UVsoft Z4 

47.2 

-8250 

U.O_o 2 

99 1 

4 1+0.02 
^•-^-0.37 

qorv + 200 

390l®°“ 

46.0l“-® 

12+11 
1 ^- 0 .3 

HE 0238-1904 Z4“ 

47.2 

-5000 

n 

^•^-0.1 

20.0l“l 

q oq+O.lO 

0.00_Q 

34001®“““ 

140l^;[ 

45.01“;® 

0 7+0-5 
' -0.2 

SDSS J0838-I-2955 

47.5 

-5000 


20.81“;® 

3.8 

3300l[®“° 

300l?™ 

45.41“;® 

9 q+1.5 

^•9_o.8 

SDSS J0318-0600 

47.7 

-4200 

-3.1 

19.9 

3.3 

6000 

120 

44.8 

0.13 


“For the high-ionization phase of trough B. 


was measured from the S iv and S iv* troughs, and one or 
both of the troughs must be saturated. Siv* cannot be sat¬ 
urated since its trough is shallower than S iv (see discussion 
in Section 6.3 of Paper I), thus the true Siv column density 
increases by a factor of eight, and the S IV*/S iv column den¬ 
sity ratio decreases by the same factor. This correction would 
decrease rie by an order of magnitude (see Figure [3|, which 
would iircrease 7? by a factor of three (see Equation ©), and 
increase M and iJk by a factor of three (see Equation ([5))). 
Although this effect would enhance the energetics of the out¬ 
flow, we suspect that S iv is not saturated since the depth of 
Siv A1062.66 is shallower than either line from the P v dou¬ 
blet which is presumed much closer to saturation. As S iv 
and Pv are both high-ionization ions, we assume that they 
have similar covering factors and thus their absorption lines 
would be the same depth if saturated. 


5.3 Comparison with other outflows 

In Table [5] we also list four other outflows with the highest 
iJk determined by our group. Two of these are low-ionization 
BALs that appear in SDSS J0838-I-2955 dMoe et al.ll2009l i 
and SDSS J0318-0600 (iDunn et al , and their distance 
estimate is based on singly ionized species. This introduces 
an uncertainty (see Section [T] for discussion) that can be 
alleviated by utilizing metastable levels from high-ionization 
ions such as Siv/Siv*. 

The first object to show a large kinetic luminosity 
(i?k=10'*® erg s“^) using Siv/Siv* absorption was SDSS 
J1106-I-1939 (Paper I). The outflow we study in this pa¬ 
per (SDSS J0831+0354) is the second such example, and it 
exhibits properties similar to SDSS J1106-I-1939, albeit at 
one-third the distance to the quasar. In a cursory exami¬ 
nation of the objects to be analyzed in our VLT/X-shooter 
survey, at least two additional BAL outflows show S iv/S iv* 
tro ughs. _ 

lArav et al.l ll2013l l also used the excited state from a 
high-ionization ion (Oiv/Oiv*) to determine the energetics 
of an outflow from quasar HE 0238-1904. Although the out¬ 
flow did not show absorption from low-ionization ions, ab¬ 
sorption from the very-high-ionization ions Nevill and Mgx 
confirmed the existence of a high-ionization phase. As dis¬ 
cussed in the previous section, the S IV and P v absorption 
in SDSS J0831-I-0354 could originate from a separate ion¬ 


ization phase, but we cannot confirm its existence without 
absorption data from additional high-ionization ions. 

5.4 Distance of quasar outflows from the central 
source 

The outflow of SDSS J0831+0354 lies closer to the cen¬ 
tral source than the other outflows listed in Table [S] How¬ 
ever, this distance (R ~ 100 pc) is three orders of mag¬ 
nitude greater than the trough forming region (0.01-0.1 
pc ) for accretion disk line-driven winds (e.g. Figu re (3) 
of [Murray et al.|[l995l : IProga. Stone fc KallmanIbOOCll) . This 
large empirical distance scale is not limited to high iJk out¬ 
flows. The large majority of measured distances to BAL 
(and narrower) outflows that were deduced using troughs 
from excited states yield distances between 10-10,000 pc 


dde Kool et alJ I 2 OOII: iHamann et al.l 

2001 

: Ide Kool et al. 

2002bja: Moe et al.l 20091: Dunn et al.l 

201 c 

: iBautista et al. 

2 OIOI: boki et al.ll201ll: lEdmonds et al. 

I 2 OIII: iBorguet et al. 

2012blal. 2 OI 3 I: Arav et al. 20131: Lucv et al. 2014|). Two out- 


flow componen ts are found to be be tween 1-10 pc from the 
central source (Ide Kool et al.|[2002bi y'). The discrepancy be¬ 
tween the empirical distances and those inferred from accre¬ 
tion disk line-driven winds models suggests that the latter 
models are not applicable to the trough formation of quasar 
outflows. 


5.5 Reliability of measurements 


Our erouD dArav 19971: lArav et all 1999bl al. l2001a bl: 

Arav. Korista & de Kooll 20021: 

lArav et al.l 20031: 

Scott et al.l 2004: Gabel et al. 

2005h and others 

(iBarlow. Hamann & Sargent 19971: 

Hamann et al.l 1997: 

Teller et al.l Il998l: IChurchill et al.l 

I 999 I: iGangulv et al.l 

I 999 I) showed that column densities derived from the 


apparent optical depth (AOD) analysis of BAL troughs are 
unreliable when non-black saturation occurs in the troughs. 
For this reason, we use partial covering (PC) absorption 
models in this analysis. 

It is instructive to assess quantitatively the difference 
in iJk that arises from using these two absorption models on 
our results. To this end we re-determine the photoionization 
solution using the AOD measurements for Pv, Sim, Aim 
and S IV (i.e. from the first column of Table[2l). The resulting 
kinetic luminosity of Eik=10'*® ® erg s“^ for the UV-Soft AZq 
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Table 6. Sensitivity of derived energetics to input changes. 


Assumption 

Effect (s) 

Outcome 

Reference 

N(Si iii)> Naod 

higher Nu, 
lower Un 

higher Ek 

Sectionl4.ll 

z < 4 Z 0 

higher Nu 

higher Ek 

Figure!^ 

different SED 

lower Nii 

lower Ek 

Figured] 

n > 0.08 

— 

higher Ek 

Equation Q 

two phases 

higher A^h, 
higher Un 

higher Ek 

Section I 5 . 2 I 

S IV* saturated 

lower riu, 
higher R 

higher Ek 

Section 15.21 


model differs only slightly from the erg s“^ that 

we derive using the PC template. Thus, our results are rather 
insensitive to the method of column density extraction. 

In Table we demonstrate that almost all possible de¬ 
viations from our representative model (UV-Soft SED, Z = 
4Zq, one ionization component, unsaturated Sim absorp¬ 
tion) results in a higher Ek for the SDSS J0831-I-0354 out¬ 
flow. Even with Ek lowered by assuming a different SED, the 
minimum energy allowed by our analysis is Ek/I/Edd = 2.6% 
for the HE0238 SED with Z = AZq (see Table [SJ. 


6 SUMMARY 

We present the analysis of BALQSO SDSS J0831-I-0354 us¬ 
ing data obtained with VLT/X-shooter. Our main findings 

are as follows: 

1. The BAL outflow in SDSS J0831-I-0354 has a velocity of 
10,800 km s“^ in the rest frame of the quasar with a width 
of 3500 km s“^ and shows kinematically connected ab¬ 
sorption from the high-ionization species Siv and Siv*. 
The ratio of the excited column density to ground was 
used to determine the number density of the outflow to 
be log (ne)=4.42±0.24 cm“®. 

2. The non-black saturation of the Pv doublet indicated 
that the absorber was inhomogenously distributed across 
the emission source. The saturation of the P v doublet 
resulted in a high ionic column density which, combined 
with the low elemental abundance of phosphorous, neces¬ 
sitates a high total hydrogen column density (Ah). 

3. Variations of the incident SED had only a small effect 
on the ionization solution ((7h,Ah) of the outflow, while 
an increased metallicty of 4Zq lowered the required hy¬ 
drogen column density by 0.4 dex. The UV-soft SED 
with Z=4Zq was adopted as the representative model, 
and the hydrogen column density was determined to be 
log (Ah)=22.0±0.4 cm“^ with an ionization parameter 
of log (Uh)=-0.25±0.3. 

4. The outflow was determined to be located 110 pc from 
the central source and contained a kinetic luminosity of 
5.2% of the Eddington luminosity of the quasar. We es¬ 
tablish half of this value (2.6%) as a conservative lower 
limit. This large kinetic luminosity makes the outflow 
of SDSS J0831-I-0354 a candidate for quasar-mode AGN 
feedback effects seen in s imulations (e.g. iHopkins et al.l 
I2OO6I : lOstriker et al.l[201ol ~). 
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